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Abstract 

We calculate the first relativistic corrections to the Kompaneets equation for the evolution of 
the photon frequency distribution brought about by Compton scattering. The Lorentz invariant 
Boltzmann equation for electron-photon scattering is first specialized to isotropic electron and 
photon distributions, the squared scattering amplitude and the energy-momentum conserving delta 
function are each expanded to order u 4 /c 4 , averages over the directions of the electron and photon 
momenta are then carried out, and finally an integration over the photon energy yields our Fokker- 
Planck equation. The Kompaneets equation, which involves only first- and second-order derivatives 
with respect to the photon energy, results from the order v 2 /c 2 terms, while the first relativistic 
corrections of order v /c introduce third- and fourth-order derivatives. We emphasize that our 
result holds when neither the electrons nor the photons are in thermal equilibrium; two effective 
temperatures characterize a general, non-thermal electron distribution. When the electrons are 
in thermal equilibrium our relativistic Fokker-Planck equation is in complete agreement with the 
most recent published results, but we both disagree with older work. 
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I. INTRODUCTION 



The Kompaneets [l| equation, 

d a T n e 1 d 4 f df{t,u) r i 1 

at^ w ) = -^7^aJ w |r- s ^- + ft[i + /(t )W )j/(t, w )| > (i.i) 

describes the scattering of unpolarized, low energy photons of frequency u on a dilute dis- 
tribution of non-relativistic electrons when all the particles — photons and electrons — are 
distributed isotropically in their momenta. The non-relativistic total photon-electron cross 
section is the Thomson cross section a T . The electron number density and mass are denoted 
by n e and m e . The photon phase space distribution f(t,u) is normalized such that the 
number n 7 of photons per unit volume is given by 

n,(t) = 2j(0±f(t,uj), (1.2) 

in which the prefactor 2 counts the number of photon polarization states and k is the photon 
wave-number vector with |k|c = u. If the electrons are in thermal equilibrium described 
by a Maxwell-Boltzmann distribution, then T is the temperature (in energy units) T e of 
this thermal distribution. However, the Kompaneets equation ( 11. ip holds for any isotropic 
distribution of electron momenta with T defined to be 2/3 of the average energy in this 
For photons with a Planck distribution, 



distribution 



a. 



^■'^M' ^W-r (L3) 



The terms in the curly braces in the Kompaneets equation (II. ip vanish when T 7 = T. In 
particular, if T = T 7 = T e , there is a time- independent photon distribution in thermal 
equilibrium with the electrons 1 . 

Our purpose here is to examine the first relativistic corrections to the Kompaneets equa- 
tion. These corrections have been previously computed by Challinor and Lasenby (C&L) 
3) for the case in which the electrons are in a thermal distribution. Using the method of 



C&L, Itoh, Kohyama, and Nozawa ^| carried out the expansion to a much higher order 



in v/c. Subsequently, Sazonov and Sunyaev [5] confirmed the previous work of Challinor 



1 Since Compton scattering preserves the photon number, the collision term on the right-hand side of 
Eq. (|l.ip also vanishes for a general Bose-Einstein distribution of massless particles at temperature T, 
/(i,w) = [exp{(^/T) - a} - l]" 1 . 
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and Lasenby. Here we use a method that is quite different from that employed by C&L, a 
method that does not require that the electrons be in thermal equilibrium. Although the 
structure of our result is quite different from that found by C&L, we agree with C&L in 
the number of higher-order derivatives with respect to the photon frequency u which must 
supplement the Kompaneets equation to correctly account for the relativistic corrections. 
Such higher-order derivative terms are missing from the ad- hoc treatments of Cooper jf}] 
and of Prasad, Shestakov, Kershaw, and Zimmerman These authors assume (incor- 
rectly) that the relativistic corrections may be accounted for by simply replacing the factor 
u A that stands just before the curly braces in Eq. ( II. ip by a function a(u),T) which is de- 
termined so as to give the rate of change of the photon energy density including the first 
relativistic corrections. We compute both the rate of energy exchange between the photons 
and electrons and the Sunyaev-Zel'dovich effect 9Mll| which follow from the relativistically 
corrected Kompaneets equation. Including the first relativistic corrections, our results en- 
tail two effective temperatures T e g- i and T e s2 which are defined by energy moments of the 
electron phase-space distribution. When the electron distribution is restricted to a thermal, 
relativistic Maxwell-Boltzmann distribution at temperature T, T e gi = T e g2 = T and we 
find, after some algebra, that our results that have a completely different structure are, in 
fact, in complete agreement with those of C&L. Moreover, the rate of energy exchange that 
we compute (also written down by C&L) agrees with that found earlier by Woodward 8]. 

Our presentation is organized as follows: After describing the general method we use in 
Sec. HIl we then outline the calculation in Sec. Hill using the results of several Appendices. 
Finally, our results are shown in Sec. HV1 Sec. IIV Al presents our general result, Sec. IIVBI gives 
its restriction to the case in which the photons are in thermal equilibrium at temperature 
Try, Sec. IIV CI derives the rate of energy transport between photons at temperature T 7 and 
the electrons in a general distribution, and finally, in Sec. IIV Dl the Sunyaev-Zel'dovich effect 
for non-thermal electrons with the first relativistic correction is briefly described. 
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II. RELATIVISTIC BOLTZMANN EQUATION FOR ISOTROPIC SCATTERING 



We start from the Lorentz invariant form of the Boltzmann equation for electron-photon 
scattering: 



Here we revert to units in which % — 1 = c, but we shall return to conventional units when 
we write the final result. The left-hand side of the equation involves the relativistic scalar 
kd = u)(d/dt) + k • V. We are assuming that the electrons and photons are not polarized. 
Hence |T| 2 denotes the square of the Lorentz invariant scattering amplitude that is summed 
over the initial and final electron and photon spins. It is divided by the initial electron spin 
weight g e = 2 so as to describe the average scattering from an initially unpolarized ensemble 
of electrons. It is divided by the square of the photon spin weight g^ = 4 because initially 
there is an unpolarized mixture and finally the scattering is into the scalar density f(x, k) 
that describes a typical photon [with the factor g> 7 = 2 needed to provide the photon number 
count in Eq. (II. 2p ], The function g(x,p) is the electron phase space density. We choose our 
Lorentz metric to have signature (— + ++) so that t = x° = —x while for the spatial 
coordinates x k — x^. 

We now specialize to the isotropic case of interest where f(x, k) — >■ fit, to) and g(x,p) — > 
g(t, E), with the electron number density given by 



The integration variables p and p' in Eq. ( 12. ip are dummy variables. We shall make the 
interchange p <H- p' in the first set of terms in Eq. ( 12.11) so as to have a common factor of 
g(t,E) for the two 'scattering in to' and 'scattering out of terms. To keep a convenient 
form, we shall also use the detailed balance relation 



for this first term in Eq. (12. ip . Finally, we note that the p' integration is best performed 




T{p', k'; p,k)\ { [1 + f{x, k)\ 2g{x, p') fix, k') - [1 + f(x, k')} 2 g{x, p) fix, k) } . 

(2.1) 




(2.2) 



\Tip',k';p,k)\ 2 = \Tip,k;p',k')\" 



(2.3) 
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using 



(d 3 p') 1 (d 3 p>) 



(dV 



(2tt)3 W (2tt)3 2v / p '2 + m 2 (2vr) 
against the four- dimensional delta function which now replaces 

p' = p + k — k' , 



5 (/ + ml) 



(2.4) 



(2.5) 



giving 



In this fashion, we obtain 



p' 2 + m 2 e = 2p(k- k') - 2kk' . 



^ f& = [ Oi Si 2 <?(*, E) 2n 5 (2p (k - k') - 2kk') 



(2tt) 3 2lo' (2tt) 3 2£ 
{|T(p',A;;p,A;')| 2 {[l + /(t,u;)]/(t,a;0 

-T(p',A;';p,A ; )| 2 [l + /(t,a;')]/(t,a;)} 



(2.6) 



(2.7) 



in which the four-momentum p' in T(p', k') is determined by Eq. (12.50 . 

The angular part of the integrations over p and k' pick out the completely rotationally 
invariant part of the integrand. Thus, with angular brackets denoting the average over all 
the orientations of the vectors within it, we may make the replacement 



S(2p (k - k') - 2kk') |T(p', k'; p, k) 



->• (5(2p (k - k') - 2kk') \T(p', k'; p, k) 
= s(p;u',u). 



(2.8) 



In view of these remarks, we may write Eq. (12.70 as 
with 



(2.9) 
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F(t,u;p) = f- / u' dJ {s(p;u,u')[l + f(t,u)]f(t,u') - s(p;u',u)[l + f(t,u')]f(t,u)}. 

Z7T JO 1 > 

(2.10) 



For the evaluation of Eq. (12.100 it is convenient to separate the weight that appears there 
into symmetric and antisymmetric parts: 



s(p; to', uj) = s s (p; uj', uj) + s A (p; uj', to) 



(2.11) 
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with 



and 



s s (p; u, u') = +s s (p; u', u) , 



s A (p; u), u') = —s A (p; u', u) . 
With this decomposition, Eq. (12.10j) becomes 

s s (p; u', u) - s A (p; u', u)] f(t, u') 



u r°° 

F(t,u;p) = — u'du' 
Ztt Jo 



f(t,u) / u'du' s s (p]u',u) + s A (p]u',u) 

2n Jo L 



u 



2f(t,u) / u'du's A (p;u',u)f(t,u'). 

2,71 Jo 



(2.12) 



(2.13) 



(2.14) 



III. EXPANSIONS AND ANGULAR AVERAGES 



It proves convenient to use the angles a , a' between p and k , k', and the angle 9 between 
k and k'. We also use the velocity 



iPl 



iPl 



E E(\p\) 
The delta function in Eq. (12.81) now becomes 



< 1 



(3.1) 



S{2p (k - k') - 2kk') = — 5(u (1 — v cos a) —u'(l — vcosa') — (uu'/E) (1 — cos#) ) 



2E 

and the squared scattering amplitude (IA8[) now appears as 

1 — cos 9 



(3.2) 



\T(p' } k'-p,k)\ z = 6rrmia T {2 



1 — v cos a) (1 — v cos a') 

2 (1-COS0) 2 



uu 



+ ( 1 - 2 ) a 



v cos a) 2 (l — v cos a') 2 



/;2 (l-cos^)-2(l- V 2 ) 

(3.3) 



As we shall see, the Kompaneets equation results from expanding Eqs. ( 13.21) and (13.31) in 
powers of v, with this equation resulting from the order v 2 terms. The leading relativistic 
corrections that concern us require that the expansion be carried out to order v 4 . We note 
that since in the applications that we envisage, u, u' ~ T ~ p 2 /m e ~ v 2 E, u/E or u'/E 
should be counted as being of order v 2 . 



6 



The needed expansion of the delta function (13. 2\i in powers of v reads 
2E 5(2p (k - k') - 2kk') = sLj - J) 

v [oj cos a — oj' cos a') 5' (ui — oj'j 

- v 2 [oj cos a — oj' cos a') 2 5" (oj — oj'\ 



+ 



'l-cos9)5' oj-oj' 



OJOJ 

~E~ 



— v 3 (oj cos a — oj' cos a') 3 5"'(oj — oj' 
— v (wcosa — a;' cos a') — — - (1 — cos^) 5"{oj — oj'^J 



+ 



— v(oj cos a — oj cos a) i — oj J 

la/ / l\2 OJOj' . . ... / ,\ 

-f (a; cos a — oj cos a ) — — (1 — cosy) o [oj — oj J 

+ 0(v 5 ) . (3.4) 



+ 1 l^r\ (1-cos6) 2 S"(oj-oj') 



2 \ E 

The term of order v°, S(oj — oj'), makes no contribution, since for it the two parts of the 
collision integral, the 'scattering in to' and the 'scattering out of, cancel. Hence we need 
expand the squared amplitude (13. 3 j) only to order v 3 to obtain results good to order v 4 : 



\T(p',k';p,k)\ 2 ~ 6nm 2 e a T \v° 



2v (1 — cos i 



+ v (1 — COS ( 



(i + cos 2 e 



(cos a + cos a') cos 9 



{cos 2 a + cos 2 a' + 2 cos a cos a' J 



+ (2 — 3 cos 2 a — 3 cos 2 a' — 4 cos a cos a'j cos # 



+ (1 - cos 



( — 2 cos a — 2 cos a' + 2 cos 3 a + 2 cos 3 a 

+4 cos 2 a cos a' + 4 cos a cos 2 a') 
— ( — 4 cos a — 4 cos a' + 4 cos 3 a + 4 cos 3 a' 

+6 cos 2 a cos a' + 6 cos a cos 2 a' ) cos 6* 



(3.5) 



We now multiply the expressions (13 .4p and (13. 5 j) together and retain the resulting terms 
up to those of order f 4 . We then use the result of angular averaging over the directions of p 



detailed in Appendix [B] and subsequently average over the direction between k and k', the 
average over cos 9. To present the results in a compact form, we separately record the order 
v 2 result, the one that gives the Kompaneets equation, 



s 2 (p;u',u) = (\T(p',k' ] p J k)\ 2 d(2p(k-k')-2kk')} 

m 2 e f 2 lolo' „ ( A 
= 4n-^a T lv 2 — 6"(u-u/) 



E 



5'{u-u') 



(3.6) 



and the order v A result which gives the leading relativistic correction to the Kompaneets 
equation 



s 4 (p;a/,w) = (\T(p',k';p, k)\ 2 d(2p(k — k') — 2kk'^) 

7T ml [ A I" , ,v2 ,14 



15 E 



e I 4 

cr T < V 



2 (a; — J wcj+— w 
o 



<T'(u;-u/) 



-10 



— f 2 (a; — a;') u> u>' + ((u> — u') 2 + —u u' 



LdLd 

~E~ 



+ 
+ 



16 a , 9 , wu' / a; a;' 

-— v A uu' + 4v 2 (cj-cj') — + 42 — - 
5 E \ E 



2 i 



S"(u-(J) 



12 



U1U1 

~E~ 



5'(u — oj'J > . 



(3.7) 



There are no terms of odd order in v in the somewhat lengthly algebra required to obtain 
these formulae. They simplify with the aid of the delta function identities presented in 
Appendix O 



2nm 2 



s|(p; a/, tu) = 2 a T t> z 



( tu tu 



3E 



(3.8) 



and 



uu 
2-nm 2 



$2 (p; u/, oj) = —2 cr a 



E 2 



A2 



w — u; 



s 4 (p; a; , a;) = cr T 4 - v uju d [u - u ) 

2tt mi 15 E 5 



+ 



5" (w - a;') 



(3.9) 



(3.10) 



2k m\ 



s£(p;a/,w) 



A2 



-cr T U 



(W) 
15 



14wa/<Hw-a/] +285' (w-u)') 



(3.11) 
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IV. RESULTS 

It is now a straightforward although tedious matter to insert the forms above into 
Eq. ( I2.14p . perform the u' integrals, and place the results into Eq. (I2.9P to secure the 
Kompaneets equation and its leading relativistic corrections. 

A. General Result 



co 2 §- t f(t,u) =a T -^-^u i { T cS x ^TZl + % [i + /(t , W )] f(t : W ) 



n„ d 



df(t,u) 



m e c du 



duj 



1 /5 2 2 \ d/(t,w) 



m e c 2 V 2 

10m e c 2 
21 rip o~ T d 2 



47 



duo 



2m P c 2 



T effl fr [l + f{t,u)\ f{t,co) 



6T eff i 



' df(t,u) Y _ h df(t,u) + t2 1 ^ 3 /(t,a; 



+ 



\ du ) du 



7 n P cr T d? 



fi dcu 3 
+ 3T cS1 h[l + f(t,u)]f(t,u) 



10m 2 c 3 do; 3 



(T 2 S2 + T 2 S1 



do; 



2T Bffl n[l + /(t,o;)]/(t,a;) 



(4.1) 



Here we have reverted to conventional units and used the effective temperature definitions 
( TD5j) and (ID7j) in Appendix |D] which, for convenience, we repeat here: 



T e ffl — 



(d 3 p) p 2 c 2 



and 



rp2 
J cfT2 



(2tt) 3 3E(p) 



2g{t,E), 



2g{t,E). 



(4.2) 



(4.3) 



15 n e J (2vr) 3 V2m e/ 

We have chosen to order the terms in our result (14.11) so as to have successive parts involve 
overall higher derivatives. Since each part starts out with at least one overall derivative, the 
result conserves photon number as it must. Only the part with the single overall derivative 
d/duj contributes to the rate of energy exchange between the photons and the electrons. 
Similarly, the rate at which the second moment (hco) 2 changes with time is affected only 
by the parts involving d/du and d 2 /du 2 while all the parts contribute to the time rate of 
change of the (hu) 3 moment. We have kept some photon frequency derivatives within the 
sequence of increasingly higher overall derivative so that the sum of the terms in each of these 



groups vanishes in thermal equilibrium. Each of these groupings in the result (14. ip vanishes, 
in fact, in the more general situation in which the photon distribution is of the Planck 
form (jl.3p but with the electron g(E) constrained only to have T c ff i = T e g 2 = T T With 
the electrons in thermal equilibrium, T e m reduces without approximation to the electron 
temperature T e , and T e g 2 — > T e in first approximation, an approximation that suffices since 
T e ff 2 only occurs in the relativistic correction terms. In this case, each of these groupings of 
terms vanishes for a Planck distribution when T 7 = T e . Our expression (14.11) is in complete 
agreement with the work of Challinor and Lasenby [3j in the limit in which the electrons are 
in thermal equilibrium at temperature T e . However, as mentioned previously, the structure 
of our expression (14. ip for electrons in thermal equilibrium — which is equivalent to that of 
C&L — differs completely from the previous (incorrect) results of Cooper jf| and of Prasad, 
Shestakov, Kershaw, and Zimmerman [?]]. 

B. Photons In Thermal Equilibrium 

For photons in thermal equilibrium at a temperature T 7 [the Planck distribution (II. 2p 
or, more generally, a Bose-Einstein distribution of photons] the result reduces to 




(4.4) 



C. Energy Transport 



The energy transfer per unit volume to the photons is given by 




(4.5) 



10 



From Eq. (I4.4p . we see that, including the first relativistic corrections, this energy transfer 
between a photon distribution in equilibrium at temperature T 7 and an arbitrary isotropic 
distribution of election energies involves 



2 J^f^ = ^i exp 



7 ) 

o | oo i o| 2 

^E^^CM^T^, (4.6) 



and 



(rf 3 k) 
(2vr) 3 



/oM = £ ~2 / exp j ^ 1 

£1 ff J ° I T 7 J 



Hence, again reverting to conventional units, 



u 7 



4 ** £ K 1 - ^ ^) ( r - - + ( T - - r -> ■ < 4 - 8 » 



In the electron thermal equilibrium limit in which T e gi = T e Q2 — T e , the energy transfer 
rate (14. 8 j) agrees with the rate given by Woodward [8] which was later confirmed by Challinor 
and Lasenby [3]. The result (14. 8 p holds, of course, for a general isotropic distribution of 
electrons with the two effective temperatures T e fn and T e fj 2 defined by the integrals (14. 21) 
and (OH. 



D. Sunyaev-Zel'dovich Effect For Non-thermal Electrons 

Equation (14. 4 p can be used to generalize published results on the Sunyaev-Zel'dovich 
effect, the distortion of the cosmic microwave background, a Planck distribution at a very 
low temperature, by high energy electrons in hot plasmas in galactic clusters. The distortion 
involves changing the time derivative in the Boltzmann equation to the proper coordinate 
distance £ along the line of sight through the plasma cloud, t —> £/c. In doing so in the 
relativistic corrected Kompaneets formula (I4.4p . we encounter two dimensionless variables 

Vl = a T [d£n e (£)^§-, (4.9) 

and 

y 2 = a T [d£n e (£) r 3^, (4.10) 
11 



that replace the conventional y parameter 2 . It is now convenient to define 

i 7 

and write the Planck distribution as fo(x) = [exp{x} — l] -1 . Since the microwave background 
temperate T 7 is so low, while the electrons have average energies that are several keV, we may 
neglect the very small ratios T 7 /T effl and T 7 /T eff2 as well as T 7 /m e c 2 . With the omission of 
these terms, carrying out the derivatives in Eq. (14. 4p yields the spectral distortion 



= Vl x [1 + f (x)} {x-A + 2 xfo(x)} 



+y 2 x [l + / (z)] (235 - 84x + 7x 2 ) - 10 



10 

+ ^ (235 - 252 x + 49 x 2 ) x f (x) 

+ 1 -f(x-2)x 2 fZ(x) + ^x 3 f (x)}. (4.12) 

In the limit of small x, the Rayleigh- Jeans region, Eq. (I4.12p simplifies to 

Af(x) 17 

In contrast to previous papers on the Sunyaev-Zel'dovich effect, Eqs. (I4.12p and (14. 13j) hold 
when the electrons are not in thermal equilibrium; when they are in equilibrium Eqs. (14. 12ft 
and (I4.13P agree with the expressions obtained by Challinor and Lasenby 3 j^|. 



2 Note that even in the case of a plasma in local thermodynamic equilibrium with an electron temperature 
T e (£), the relativistic corrections involve a different y parameter (y = y 2 ) defined by the electron number 
weighted average of T^(£) rather than the first power T e (£) that appears in y — y\. 

3 These authors, however, do not define the proper parameters y\ and yi that, as noted in the previous 
footnote, are needed for the relativistic treatment, but rather use a y = y\ and then multiply this by an 
undefined electron temperature to obtain a yi parameter for the relativistic corrections. 
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Appendix A: Squared Amplitude Details 



We first express the fully relativistic squared amplitude [12| as 



\T(p', k';p, k)\ 2 = Qitml 



K K \ 

- + - +2 

K K I 



' m 2 e m 2 e 

K K 1 







( ml 




{- 


—7) j 



(Al) 



in which a T = 87rrQ/3 is the Thomson cross section with r the classical electron radius and, 
with our space-like metric 



k = —pk = p°k° — p • k , k' — —pk' = p°k'° — p • k' . 

Note that in terms of these variables the delta function constraint ( 13. 2 j) reads 

k' — k = kk' . 

It is convenient to use the variable 



k = y KK' 



so that the relation (IA3j) may be written as 



k' kk' 
— = 1 + — 

K K 



! k' kk' 

K K 



The proper solution of this quadratic equation, written in terms of k'/k, is 

kk' 



Similarly, 



k!_ _ 1 
~k ~ 2 



K 1 

k/ ~ 2 



klf 

K 



+ 2 + 



(kk r 



+ 4 



K 



'kk'Y kk' 

v « ) R \ 



x) +4 



Making use of Eq. ( 1A3I) and some algebra now presents 



kk' 



'kk'" 



\T(p\k'-p } k)\ z = Qnm 2 e a T {[ — ) +2+ 2m^- + ( — 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



Appendix B: Angular Averages 

The calculation outlined in the text involves the integration of 

cosa = p ■ k = p l k l , cos a' = p ■ k' = p l k a , 



(Bl) 
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and of the products of the powers cos m a cos n a', over the solid angle of the momentum p 
associated with the electron distribution g(t,E). This is equivalent to averaging over all 
orientations of the unit vector p, and this averaging can be performed at any stage of the 
computation — it may be performed before the actual integral over p is carried out. The 
averages may be expressed as contractions of outer products k l ■ ■ ■ and k' m ■ ■ ■ with the 
rotationally invariant tensors that result from the angular averages {p k p l • • •). For example, 

(cos a cos a') = U k' m (p k p m ) . (B2) 

This is the method that we shall employ. 
Under the angular average 

(p l ) = 0, (B3) 
and since p is a vector, not a pseudo-vector, 

(p l p m p n ) = 0, (B4) 

because the only rotationally invariant, third rank tensor, is the pseudo-tensor e lmn . The 
lowest-order correlation is 

<p l p m ) = ±SF», (B5) 

where 5 kl is the matrix element of the invariant unit matrix. The overall coefficient is 
determined by the trace (p l p l ) = (1) = 1 . The final average that we shall need is 

(p k p l p m p n ) = [d kl 5 mn + 5 km 5 ln + 5 kn 5 lm ] . (B6) 

Here the particular combination of the delta symbols on the right-hand side is required to 
reproduce the complete symmetry of the left-hand side under any permutation of the indices 
k , I , m , n. The overall coefficient again follows from taking the trace over any index pair 
and comparing the result with Eq. (IB5j) . 
Therefore, 

(cos a) = = (cos a') , (B7) 



and 



from which follow 



(cos a cos a) = - k • k' = - cos 8 , (B8) 



(cos 2 a) = \ = (cos 2 a') . (B9) 
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Since the angular average of three momentum vectors vanishes, 

= (cos 3 a) = (cos 2 a cos a') = (cos a cos 2 a') = (cos 3 a) 



(BIO) 



Next, 



(cos 3 a cos a') 
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k 2 kk' + k 2 kk' + kk'k 2 



and 



cos a cos 3 a') = - cos 9 . 



cos#, (Bll) 



(B12) 



Finally 



15 



cos 2 acos 2 c/> = — k 2 k' 2 + (k-k') 2 + (k-k') 2 =— l + 2cos 2 # , (B13) 



15 



from which follow 



(cos 4 a) = - = (cos 4 a'} . 
5 



(B14) 



Appendix C: Delta Function Identities 

The work in the text involves various derivatives of 6(u—u') = 5(x) multiplied by various 
powers of u — u' = x. Simple manipulations can be performed to place a derivative d n /dx n 
to the left of a power x m , leaving lower derivatives and lower powers of x. Again, the lower 
derivatives can be ordered to the left so that all the derivatives appear as total derivatives. 
According to the rules of generalized functions, any resulting term of the form x l 8(x) gives 
a vanishing contribution. The simplest example of this procedure is 

x5'(x) = 4- \x6(x)]-S(x) . (CI) 
ax 

As we have just noted, the first term may be discarded. Moreover, the delta function 
S(x) = S {uj — uj') with no derivative can also be omitted because it gives rise to equal 
contributions from the 'scattering in to' and 'scattering out of terms in the Boltzmann 
equation which cancel. We use the symbol ' =' to denote the only effective parts that 
remain after the manipulations described above have been made. Thus, we write Eq. (ICip 
as 

x5\x) l =' 0. (C2) 
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More involved computations lead to the effective results 

x8"{x) '=' -2 8'(x), 



(C3) 



x 2 S"(x) '=' , 



(C4) 



x8"\x) 



3 8"(x), 



(C5) 



x 2 8'"(x) '=' 6<5'(x), 



(C6) 



x 2 8""{x) '=' +12<S"(a;), 



(C7) 



and 



x 4 8""(x) '=' 0. 



(C8) 



Appendix D: Effective Temperatures Denned By Electron Distribution Integrals 

Here we shall explain the definitions of the two effective temperatures T e gi and T e gf 2 that 
reduce to the electron temperature T e when the electron relativistic phase-space distribution 
g(t,E) is restricted to be a Maxwell-Boltzmann distribution. 

For our system of free, relativistic electrons, the number density is given for an arbitrary 
phase space distribution g(t,E) by 




(Dl) 



For the case of thermal equilibrium, 




(D2) 



in which 



E(p) = p 2 + m 2 e 



21 V2 



(D3) 



is the total relativistic energy of an electron with momentum p. 
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To obtain the definition of T eff i, we note that partial integration gives 4 



n, 



(d 3 p) 



(2tt) 3 
^ (d'p) 



2 const, exp 



T J (2n 



E(P) \ ld_ 

3 dp P 

r E( P )\ 



T 



— const, ip- A 2 exp r — ^ 



(d 3 p) p 2 



2 </(*,£). 



Thus we define 



cff 



(2vr) 3 3£(p) 

(d 3 p) p 2 



- 1 / 



2g(t,E). 



(D4) 



(D5) 



(2tt) 3 3£(p) 

We have just shown that in thermal equilibrium, T cffl — > T is an exact, relativistic result. 

To obtain the definition of T eff2 , we note that we are computing only the first relativistic 
corrections to the Compton Fokker-Planck equation. Hence, for an equilibrium distribution, 
we may approximate a relativistic correction integral by 

2 



2,(^)^2 cons, exp { -. e /T}/g 



(rf 3 p) 

) 3 



P 



2m P 



exp 



P 



2m e T 
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Therefore, for an arbitrary distribution, we shall define 

j3 * ^ / 9 \ 2 



(D6) 



'cff 2 



(d* P ) 



P 



15 n e y (2tt) 3 V2m e 



2 



(D7) 



To leading order, T cS2 — > T when i?) becomes a thermal distribution. 

With these results in hand, we can now evaluate the the integrals needed in the text. 
The following two integrals appear with the lowest-order functions s| and s£ , and thus they 
must be evaluated to both lowest and first non-leading orders: 



(d 3 p) v 2 
(2tt) 3 E 2 



(d 3 p) p 2 
(2tt) 3 E 4 

(d 3 p) 1 



g(t,E) 



(27r) 3 m 
3 n f 



p 
~E 



2 m\ 



T e s\ — 



6 

m e 

15 T 2 ff2 
2 m P 



p 



2m, 



(D8) 



4 What follows is a proof that, for a dilute, non-interacting gas, the familiar equation of state p = nT holds 
even for relativistic particles. 
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(d 3 p) 1 



J (2tc 



(2tt) 3 E 3 



g(t,E) 



f(d 3 p) 1 




/ (2tt) 3 ml 


2 m\ 



g(t,E) 



2m 3 



9 T effl 



2 m P 



(D9) 



On the other hand, for the higher order s| and S4 terms, we need only leading evaluations: 



(2tt) 3 £ 2 y{ ' ; ~ 2 m 2 ' 



y f27T 



^ — oft 



2mj 



(d 3 p) v 2 3 n e 



(DIO) 

(Dll) 
(D12) 
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